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Abstract 



We propose that extra dimensions might be detected with higher dimensional 
analogues of black holes. The usual 4-dimensional acoustic (sonic) black hole metric is 
extended to arbitrary dimensions. The absorption cross-section of Hawking radiation 
\ on the brane and in the bulk are calculated in the semiclassical approximation. 
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1 Introduction 



A lot of interest in recent years has been raised for field theories where the 
standard model of high-energy physics is assumed to live on a 3-brane embed- 
^ ' ded in a larger space-time, while only the gravitational fields are in contrast 

usually considered to live in the whole spacetime [TII^II?|H] . Arkani-Hamed et al 
[2], proposed a new framework for solving the hierarchy problem which does 
not rely on supersymmetry or techi-color. The novelty in this idea was that 
the traditional picture of Planck-length-sized additional spacelike dimensions 
{Ip ~ 10~^^cm) was abandoned, and the extra dimensions could have a size 
as large as 1 mm. The fact that we do not see experimental signs of the ex- 
tra dimensions despite that the compactification scale of the extra dimensions 
/ic ~ l/Ki^''" would have to be much smaller than the weak scale, implies that 
only gravity can propagate in the extra-dimensional spacetime and all ordi- 
nary matter: electromagnetic, weak and strong forces, is restricted to live on 
a (3+1) dimensional hypersurfaces, a 3-brane. 

The presence of extra dimensions in brane world gravity models will 
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inevitably change the properties and physics of black holes. In brane-world 
scenario, the production cross-section for black holes is greatly enhanced. De- 
tectable signals of black holes might be found via Hawking evaporation to 
brane-localized modes. If the scale of quantum gravity is near a Tev, the 
Large Hadron Collider (LHC) will become a "black hole factory" [5] . The cor- 
relations between the black hole mass and its temperature, deduced from the 
energy spectrum of the decay products, can test Hawking's evaporation law. 
The absorption cross-section (and also grey body factors) of black holes in dif- 
ferent physical scenario have been investigated by several authors piTfSI^ITU] . 
The emission rate of 4-dimensional acoustic black holes for S wave has been 
studied in Ref. [11] . In this paper, we would like to investigate the absorption 
cross-section of higher dimensional acoustic black holes. 

The remarkable work of Unruh, in 1981, developed a way of mapping cer- 
tain aspects of black holes in supersonic flows and pointed out that propagation 
of sound in a fluid or gas turning supersonic [12] , is similar to the propagation 
of a scalar field close to a black hole, and thus experimental investigation of 
the Hawking radiation is possible. From then on, several candidates have been 
considered for the experimental test of the analogue of black holes |13j . 

In the present study, we extend 4-dimensional acoustic black hole metric 
into arbitrary dimensional space-time in Section 2. Then following the meth- 
ods of Ref. [9], we calculate the corresponding absorption probabilities and 
energy emission rates of such black hole analogues. In section 3, we present 
the calculation of bulk emission in (4 -t- n)-dimensions for / > case by using 
the low energy perturbation method. Section 4 performs the calculation of 
brane-localized scalar emission, where the brane is embedded in a (4 + n)- 
dimensional bulk. We give the conclusions in Section 5. 



2 Metric of higher dimensional black hole analogues 

Since the dynamics of effective black hole analogues metric are not de- 
scribed by the laws of gravity (i.e. Einstein equations) in general, to extend 
the four dimensional effective metric we should start with the equations of 
fluid dynamics. We consider one kind of fluids which is permitted to propa- 
gate in (4 + n)-dimensional space-time. The fluid might be condensate states 
of gravitons or some higher- dimensional scalar particles. In the following of 
our discussions, we simply neglect quantum effects of such special quantum 
fluids and assume they obey the classical fluids equation (Euler equation) in 
the low energy approximation. The fundamental equations of fluid dynamics 
in (4 -|- n)-dimensional flat spacetime are the equation of continuity 



'iP + V ■ ipv) = 0, 
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and Euler's equations 

dv -» 
P^ = p[dtv + {v-V)v]=F, (2) 

where V = + ey-§^ + + ... + e„^, and F = -Vp - pV0. Here 
denotes the Newtonian gravitational potentiaL It is also assumed that the fluid 
to be irrotational and inviscid, which imply V xv = 0. These equations can be 
linearized in the vicinity of some mean flow solution with p = pQ + epi + Ole"^), 
P = Po + ^Pi + 0{e^),(p = (po + Eipi + 0{e'^), redefining the fields as Vh = 
and V = V<f. One can finally obtain the wave equation [H], 



dp ^ ^ , 

g^poidtipi + Vo ■ \/Lp) 



+ V 



PoVifi - ^povoidtifii + vo ■ Vifi] 



=(8) 



The above equation is identified with a massless scalar field equation describing 
the sound wave in the curved spacetime background 



1 



(4) 



with the background metric, g^^^, ~ (^) " 



2 

+2 



which is a (n + 



4) X (n + 4) matrix, where the local speed of sound is defined by = and 
i,j = l...n. In spherical coordinates, assuming Vr 7^ 0, ve^ = ve^ = .. 
then have 



0, we 



P0\ "+2 



-c'{l - %)dT' + (1 - %r'dr^ + r^nl^^ 



(5) 



which is similar to the n-dimensional Schwarzschild black hole metric [15j, 



-c 



H 



71+1 



dt^ +\l 



r 



n+l 
H 



2„ri+l 



dr' 



(6) 



where C is the light velocity and is the event horizon radius. The properties 
of higher dimensional Schwarzschild black holes in brane- world scenario have 
been discussed in Ref. |16]. However, the properties of acoustic black holes do 
not necessarily relate to the outside geometry of space-time since the dynamics 
of acoustic black holes do not obey Einstein equation. In Ref. , we have dis- 
cussed black hole analogues in brane-world scenario under some assumptions 
and found the properties of acoustic black holes in brane-world scenario are 
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similar to real black holes. If p is time and position independent, the continu- 
ity equation V ■ = then implies that Vr oc -kpi- Because of the barotropic 
assumption, p is position independent implying the pressure p and the speed 
of sound c are also position independent. We can define a normalization con- 
stant To and set Vr — C/r""*"^ with C — ctq'^'^. tq is a parameter which can be 
determined by experiments. The metric can be rewritten as 



/ ^2n+i\ ( 2n+4\ 



(7) 



2 

where ds^ = (^)"+2(is^. The temperature of an acoustic(sonic) black hole is 



dr 



given by T = ^ , where Vr = c marks the exact position of the event 



Vr=C 



horizon. For 4-dimensional cases, Vr — c^, we have the numerical expression, 
T = 6 X IQ-^Kfc/SOOm/sec] [1 mm/ro]. (8) 

For (4+n)-dimensional cases, Vr — c^:^, the temperature can then be written 



as 



T = 3(n + 2) X 10-^is:[c/300m/sec][lmm/ro]. (9) 

It is clear from the above equations that unless c is very large (i.e. to 
be the velocity of light c = 3 x 10^ m/sec), the experimental verification of 
above acoustic Hawking temperature will be rather difficult. However, the 
higher dimensional black hole analogue metric presents us an alternative way 
to detect extra dimensions other than in LHC. By detecting signals via the 
acoustic black hole's evaporation to brane-localized modes, one can in principle 
determine the exact dimensions of space-time. However, since the Standard 
Model particles must to be located to an ordinary 4-dimensional spacetime, 
extra dimensions can only be probed through the gravitational force. The 
continuity equation V • = in higher dimensions should only apply to 
gravitational waves or some unknown scalar particles. 



3 Bulk scalar emission: S wave and I > 



In this section and the next section, we shall calculate the decay rate of 
higher dimensional acoustic black holes from the near-horizon low energy dy- 
namics, which is expected to give an experimental suggestion of detecting the 
thermal radiation. 

The radiation of higher dimensional acoustic black holes is usually de- 
scribed as thermal spectrum in character with a temperature T4^n- The energy 
emitted per unit time (power spectrum) by gravity-wave black hole analogues 
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for a higher number of dimensions can be given by 

~ir - V^''"^"'^ea;p(o;/T4+„)-l(27r)"+3' ^^^^ 

where / is the angular momentum quantum number and \ k \= uj for mass- 
less particles. However, considering the nontrivial metric in the region exterior 
to the horizon, there exists an effective potential barrier in this exterior re- 
gion. This potential barrier backscatters a part of the outgoing radiation back 
into the black hole. Thus the original blackbody radiation is modified by a 
frequency dependent filtering function ai{uj), caused by the gravitational po- 
tential of the black hole, which is called the "greybody factor". Greybody 
factors are important theoretically and experimentally in that they depend on 
the number of extra dimensions and encode information on the near-horizon 
structure of black holes and can be used to identify a black hole event. 

For S wave bulk scalar emission of the (4 -|- n) -dimensional gravity- wave 
black hole analogues, the calculation is similar to that of higher dimensional 
Schwarzschild cases. The calculations do not depend on the concrete form 
of the metric. In fact, it has been proved that all spherically symmetric black 
holes, regardless of the theory in which they arise, the low energy cross section 
for massless minimally coupled scalars is always the area of the horizon |18j . 
Following the method of Ref. [10], we can obtain the spherical wave absorption 
probability, 

2Ko)»+^.m[7r(n + l)/2]r(if^) 
I I - ^^(^^^T^ f(3±^ (11) 

The greybody factor can be computed by first evaluating the absorption prob- 
ability, I A{uj) p, from the ratio of the in-going flux at the future horizon to the 
incoming flux from past infinity with boundary condition that there is no out- 
going flux at the horizon, and then using the generalized (4 + ?T,)-dimensional 
optical theorem relation 



between the absorption cross section (Ti{lj) and the absorption probability 
I A{u!) 1^ for the l-th partial waves. 

We will then derive the scalar decay modes which are not spherically 
symmetric, / ^ 0. We start with the 4 -|- n-dimensional black hole analogue by 
rewriting Eq.([7]) in the following form, 

ds'^ = -h{r)dt^ + h{r)-^dP + r'^d^l^^, (13) 

where h[r) = 1 — ( y j and we have assumed po = c = 1. The scalar wave 
equation in this background is separable if we make the ansatz (j)(t, r, 9, ip) = 
e~'"^^ Ru,i{r)Yi{fl) , where Yi{fl) are now the (3-fn)-spatial dimensional spherical 
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harmonic functions. We can then obtain the radial part equation by substi- 
tuting the above ansatz into the Eq.(jl]), which reads, 



h{r) d 

y,n+2 



h{r)r 



n+2 



dR 

dr 



+ 



h{r)l{l + n + 1) 



i? = 0. 



(14) 



The idea of [7] is to solve this equation approximately in three regions: near- 
horizon regions, far field regions and intermediate regions, and match the 
solution across the boundaries of the regions. We will use the same method 
here, but simply solve the radial equation in near-horizon and far field regions, 
and assume that keep only the lowest order terms in uo in each region is enough. 
From the change of variables 0<r— i>/i<l, we can write the scalar field 
equation Eq. (fT4|) in the form. 



,,d^R 



dK^ 



3n + 7. 
2n + 4 



dR 
dh 



+ 



2^2 



_ _ l{l + n + l) 
(2n + 4)2/i(l - /i) (2n + 4)2(1 - /i) 







i?= 0, 



(15) 



where we set {ujrY to be (curo)^ near the horizon. By redefining R{h) = h°'{l- 



h)^F{h) and removing singularities at h = and h 
can be reduced to a hypergeometric equation with a - 
and c = 2q; -|- 1, where. 



1, the above equation 



a + p + 



n+3 L 
2n+4' 



a + P 



Thus, we have, 

(fjp dp 
h{l-h)-^ + [c-{l + a + h)h] — -ahF = 0, (17) 

which has a solution the hypergeometric function -F(a, 6, c; h) [20]. The crite- 
rion for the convergence of the hypergeometric function demands that Re{c — 
a — b) > 0, which force us to choose j3 = j3-. Then, the general solution of 
Eq.(??) is, 

RNH{h) = A.h'^^il - hf F{a,b, c; h) 

+ A+/i±"± (1 - hfF{a - c+l,6-c+l,2-c;/i) (18) 

Expanding the above solution in the near-horizon region in the limit r — tq, 
or /i — > 0, and choosing a = a_, we obtain the result. 
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,."+2, 



(19) 



where y is defined hy y = - To calculate the greybody factor, we must 

impose the boundary condition that near the horizon the solution is purely 
ingoing and then we set ^4+ = [S]. 

By using the h {1 — h) transformation of hypergeo metric functions, 



F(a, 6, c; h) = ^ F{a, b; a + b - c + 1;1 - h) 

1 (c — a)L [c — b) 

+ (1 - /,)^-^-^ ^(")^(^+^-") f(c -a;c~b,c-a-b + l;l-h), (20) 
T{a)T{b) 

the near-horizon solution f|T8l) expanded in terms of 1 — is given by. 



[1 - hy 



r(l + 2a)r(l -2/3- ^) 



r(l + a-/5-ft^)r(l + a-/3) 



_ ,.-.-^ r(i + 2.)r(2/? + ia-i) 

^ ^ r(a + /3 + ft|)r(a + /3) 



F(c — a,c — 6, c — a — 6+l;l — /i) 



(21) 



where has been set to vanish. Expanding the above expression in the limit 
/i — > 1, we get. 



RnhW 



r(l + 2a)r(l-2/3 



n+3 ^ 
2n+4 > 



n+3 



,.^oy+n+i r(l + 2«)r(2/j+,„^^ 
^r' r(a + /3+|±|)r(a + /3) 



(22) 



The derivation of the far-field zone which is defined by r ^ vq. In this 
limit, h{r) ~ 1 and, by setting R{r) = /(r)/r^""^-^^/^. Eg. (1141) can be rewritten 

as 



1 df 
— H + 

dr"^ r dr 



{2l + n + l] 



/ = 0, 



(23) 



which is a (/+^^)-th order Bessel equation. The solution of the above equation 
are the Bessel functions 



B B 



(24) 
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where J;_|_n+i (cot) and Yi_^_n+i[ujr) are the Bessel functions of the first and 
second kind, respectively. For small cur <^ 1, we can expand the above formula 
into, 



B+r' B_ /2\'+^r(/ + ^) , , 



In order to match the solutions across the boundaries of near-horizon field and 
far-field zone, we need to rewrite the near-horizon solution in terms of [l — h), 
before expanding the solution in the limit r ^ vq. Matching the two solutions 
Eqs.(l22l) and ( l25l) . we obtain the ratio. 



^ / 2 \2i+n+l 

T{1 + n±l)2(/ + n±l)T(i - 2/3 - ^jTja + /?)r(« + /3 + 

w+3 I "+3 



7rr(l + a - /3)r(l + a - /3 - |±^|)r(2/5 + ^ - 1 



(26) 



The ingoing and outgoing waves of Eq. fl25|) can be decomposed by introducing 
the redefinition of amplitudes. We expand Eq. (p5|) in the limit r oo, and 
then we obtain. 
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Figure 2. Analytical results for the energy rates for scalars from a 
(4+n)-dimensional black hole in the bulk for Z = 1, where x = ur^ and y = [r^^] . 



d(oo) _ 4(00) 

1 L -^'i'n 



-t- ^, 



^out 



yjr^n + 2) sjr^n + 2) 

where A^^^^ and A^^} is defined by, 



(27) 



^(00) ^ -B+ + 

B: - iB_ 



(00) 

out 



gi7r(Z+f+l)/2 



,-m(«+2+l)/2 



V2 



(28) 



The reflection coefficient TZ is defined as the ratio of the outgoing amphtude 
over the incoming amplitude at infinity, 



n 



A 



(00) 



out 



(00) 



(29) 



The absorption probability can be written, in terms of B — B+/B^, as 

2i{B* - B) 



A{uj) 



BB* + i{B* -B) + l' 



(30) 



In the limit of wtq -C 1 and BB* ^ i[B* — B) 3> 1, the absorption probability 
can approximately be written as 
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, ,2 87r(2n + 5-2/) 



(/ + B±i)(2/-1)2 



^ ^^0 \2;+n+2 '^4 2n+4/ U(n+2) 2n+4>' x 

. 2^ r(/ + H±img-^)2 



From Figure 1, we can see that if we fix the angular momentum number 
and vary only the number of extra dimensions, the absorption probability 
decreases as n increases, since the expansions of | ^ ^ is in powers of uvq <^ 1. 
Thus, I ^ P should become more and more suppressed as n increases. The 
same behavior is observed if we fix instead n and vary /. 

According to Figure 2, the emission rate of scalar fields in the bulk is 
enhanced as the number of extra dimensions increases. This is caused by the 
increase of the temperature of gravity-wave black hole analogues, which finally 
overcome the decreases in the value of the greybody factor and causes the 
enhancement of the emission rate with n at high energies. 



4 Brane-localized scalar emission for / > 



If the acoustic black hole is formed from matter on the brane, the scalar 
field is confined on the 3-brane embedded in a (4 + n)-dimensional space-time. 
The induced metric on the brane will be, 



h{r)dt' + h{r)-^dr' + r\de' + sin' 



(32) 



where h{r) = 1 — (ni)2{"+2)_ q^i the brane, the event horizon is still at r = tq and 
its area is = ATrr^. This induced acoustic metric on the brane is certainly 
not the 4-dimensional acoustic geometry. The calculation of Hawking radiation 
relies on mainly on properties of the horizon, such as its surface gravity. We 
shall calculate the absorption cross-section of the brane-localized scalar field, 
which according to Ref. [21j, most of the energy radiated by black holes goes 
into modes on the brane. Using the separation of variables, 0(t, r, 6, ip) = 
e~^^^ Rwi{r)Yi{fl) , where Yj(fi) are now the usual three-dimensional spherical 
harmonic functions, the radial equation of Eq. (jlj) is written as. 



h{r) d 

J.2 f^j. 



h{r)r' 



,dR 
dr 



+ 



= 



(33) 



Similar to the discussions in section 3, we solve this equation in two regions: 
the near-horizon region and far-field region. Our starting point is the Klein- 
Gordon equation under the near-horizon metric background. In terms of h, 
the radial differential equation now takes the form. 
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, , , , (PR 
h{l-h)-— + 



An + 7 
2{n + 2) 



+ 



dR 
dh 
1(1 + 1) 



(2n + 4)2/i(l - /i) (2n + 4)2(1 - /i) 



i?= 



(34) 



If we further define R{h) = h°'{l — h)^F{h) to remove singularities at points 
h = and h = 1, the above equation assumes the standard form of a hyper- 
geometric equation 



/i(l 



dh 



abF = 0, 



with indices a = a + P + j^fi = a + P and c = 1 + 2a, where a± 



2?i+4' 



(35) 



2n+2 



and (3 = (1 ± y (2/ + 1)^ — Au'^r'^). The criterion for the convergence of 

the horizon and imposing the boundary condition that only incoming waves 
exist near r ~ tq, one can find that ^4+ = for a = a_. To express the form 
of the solution for small h, we express the F{h) in terms of 1 — /i, by using the 
hypergeometric relation Eq. (!20|) . Thus after expanding R^nih) for r ^ vq, 
we find that the desired solution for near-horizon region is 



RnhW ^ 



r(l + 2a)r(l -2/3 



2n+3 ' 



r(i 



a 



P 



2n+3 
271+4 



)T{l + a + p) 



,ro,,+ir(l + 2a)r(2/3+^-ll 



2n+4 



r(« - + iTl)r(« + 



For the far-field region which is defined by r ^ r^, noting that h{r) 
setting R{r) = f{r)/r^/'^, we rewrite Eq.(1551) as, 



(36) 
1 and 



cPf 1 df 
— - H + 

dr^ r dr 



{21 + 1)- 

4r2 



/ = o. 



(37) 



The far field region solution can be expressed in terms of the Bessel function 
Jz+i/2(a;r) and Neumann function Y;+i/2(a;r), 



B. 



(38) 



At this stage, we expand the general solution Eq. fl38|) . in the low energy limit 
ujr and find that. 



R 



B.r^ 



FF 



r(/ + 3/2) V2 



UJ 



1+1/2 



Q /2y+i/2 r(/ + 1/2) 



TT 



(39) 
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Matching the solution Eq. (l36|) with (139|) . we obtain the ratio, 



21+1 



B_ \ujrQ, 
T{1 + |)^(/ + i)r(l -2P- t±3)r(a + p)T{a + /3 + igf ) 
7rr(l + a - /3)r(l + a - /5 - g|)r(2/3 + i±| - 1) 



(40) 



We can now compute the absorption probabihty | A{uj) \ ^ in the low energy 
limit utq -C 1, which goes as, 



' ^ ^ ' (2/ + i)H 2 ; r(/ + i)2r(i + ^^'^ 

Figure 3 demonstrates that if we keep n fixed and varying the absorption 
probability decreases as / increases. The dominant term becomes more and 
more suppressed by extra powers of wr^ and its numerical coefficients also 
decreases. If we fix instead / and vary n, a different behavior from the one 
observed in the case with a bulk scalar field emerges, that is to say the leading 
term remains the same since it is n-independent. Figure 4 is to compare the 
results of the absorption probability derived in the case of a brane-localized 
{n > 0) scalar field with a purely 4-dimensional {n = 0) scalar field. For higher 
partial waves (/ > 0), the value of the absorption probability in the case of 
a brane-localized {n > 0) scalar field is larger than the one for a purely 4- 
dimensional {n = 0) field, while for / = they are the same. Figure 5 depicts 
the behavior of the energy emission rates for particles with the angular number 
/ = 1 in the low and intermediate energy regime. The figure shows that the 
energy and the number of particles, emitted per unit time and energy interval 
is strongly enhanced, as n increases since the temperature of the black hole 
is given by the relation T^+n oc (n + 2)/27rro, which indicates that for fixed 
To, the temperature of the gravity- wave black hole analogues increases as n 
increases. This means the energy of the emission particles also increases. 



5 Conclusions 



In summary, we have extended the 4-dimensional acoustic black hole met- 
ric to higher dimensions that is similar to higher dimensional Schwarzschild 
metric in form but is not exactly any of the standard geometries typically con- 
sidered in general relativity. The fiuids here have been assumed to fill all the 
spacial dimensions including extra dimensions. We emphasize that the higher 
dimensional acoustic black holes discussed above is just a model, which may 
help us understand the physics beyond the standard model and are falsifiable 
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Figure 3. Analytical results for the absorption probability for a (4 + n) black holes 
analogues on the brane for n = 3, where x = ur^ and y =| .4 p. 




Figure 4. Analytical results for the absorption probability for a (4 + n) black holes 
analogues on the brane for I = 3, where x = lotq and y =\ A\^. 

by the future experiments. Although it might be difficult to detect the Hawk- 
ing effects of these analogues in experiments since the Hawking temperature in 
fact is very low, they provide us an otherwise way to probe extra dimensions. 
The scalar emission of Hawking particles in both (4 + n)-dimensional bulk 



13 



0.4 n 






1=1 


n=4 / 




0.3 










/ 

/ 


y 0.2 






/ 


n=3 




0.1 






/ 


y'' 

n=2 




00.2 


0.3 


0.4 0.5 


0.6 0.7 


0.8 








X 







Figure 5. Analytical results for the energy rates for scalars from a 
(4 + n)-dimensional black hole on the brane for Z = 1, where x = lovq and 

scalar field and 4- dimensional brane-localized scalar field were studied respec- 
tively in a higher dimensional gravity-wave black hole analogues background. 
The amplitude probability in a bulk scalar field was obtained that allows one 
to find its dependence on the number of extra dimensions n and the angular 
momentum number I. We found that if we fix the angular momentum number 
and vary only the number of extra dimensions, the absorption probability de- 
creases as n increases, and | ^ P should become more and more suppressed as 
n increases. The same behavior is observed if we fix instead n and vary /. The 
case in which the scalar field is confined on a 3-brane in a higher dimensional 
spacetime background is also discussed. The resulting absorption probability 
depends only on the angular momentum number through (a;ro)^'"'"^. If we keep 
n fixed and varying /, the absorption probability decreases as / increases. But, 
in both bulk and brane-localized cases, the energy emission rates are enhanced 
as n increases since the temperature of the gravity- wave black hole analogues 
increases as n increases. 
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